Optimal Surface-Electrode Trap Lattices for Quantum Simulation with Trapped Ions 
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Trapped ions offer long internal state (spin) coherence times and strong inter-particle interactions 
mediated by the Coulomb force. This makes them interesting candidates for quantum simulation 
of coupled lattices. To this end it is desirable to be able to trap ions in arbitrary conformations 
with precisely controlled local potentials. We provide a general method for optimizing periodic 
planar radio-frequency electrodes for generating ion trapping potentials with specified trap locations 
and curvatures above the electrode plane. A linear-programming algorithm guarantees globally 
optimal electrode shapes that require only a single radio-frequency voltage source for operation. The 
optimization method produces final electrode shapes that are smooth and exhibit low fragmentation. 
Such characteristics are desirable for practical fabrication of surface-electrode trap lattices. 
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Trapped ions are a promising system for quantum com- 
putation and quantum simulation [Tl [2]. For the lat- 
ter, simulations of coupled lattices that are ubiquitous in 
condensed-matter systems are of particular interest, as a 
"proving ground" for establishing the viability of large- 
scale quantum simulations as well as for studying poorly 
understood systems in great detail. 

Trapping of ions over micro-fabricated planar surface- 
electrode structures has been demonstrated [3 , and the 
lithographic manufacturing process is well suited for pro- 
ducing complex electrode arrays 0^ . Further, proposed 
schemes for local control over ion interactions using mag- 
netic fields [3 16] instead of laser beams have the poten- 
tial to greatly simplify experimental setups. However, a 
dense lattice of microtraps cannot be constructed by just 
placing individual microtrap electrode sets side by side, 
because their electric fields overlap significantly and dis- 
tort their microtraps. A regular dense lattice of ions for a 
two-dimensional condensed-matter simulation should be 
located at a height z as far from the electrode plane as 
possible, because heating effects induced by the electrode 
surfaces scale strongly with this distance (proportional to 
for Johnson noise |7] and to z^^ for anomalous heat- 
ing [8]). Moreover, this distance is crucial for directing 
laser beams onto the ions to induce spin-dependent in- 
teractions. The aforementioned cross-talk effects become 
prominent as the ion-surface distance z is increased to 
about half the inter-ion spacing d. In this regime the 
surface electrodes for individual traps must be replaced 
by complex electrode patterns to generate a desired mi- 
crotrap lattice. We present an algorithm that directly 
produces the periodic electrode structure for a desired 
lattice of trapping sites. 

We consider ions of mass m and charge q that are con- 
fined by a time-averaged ponderomotive potential cre- 
ated by an inhomogeneous radio-frequency (rf) electric 
field with amplitude E{r) = — V<i>(r^ at angular fre- 



quency riff. Additional electrodes inducing static electric 
fields could be incorporated for greater fiexibility of de- 
sign, but we do not consider this possibility here. In the 
adiabatic approximation, which assumes that the motion 
of the ions is slow on the time scale 2Tr/Q,^f, the pondero- 
motive pseudopotential is [3 EHl E] 



(1) 



Our aim is to create a periodic pseudopotential lattice 
with minima at a set of desired trap positions rj. We 
would like to specify the curvature of the pseudopoten- 
tial at these sites, and also require that the rf field E{fj) 
vanishes, so that the ions undergo no rf-induced micro- 
motion [10 . At a field-free point rj, the pseudopotential 
curvature tensor is proportional to the square of the elec- 
tric potential curvature tensor ^'^•'(r*,) — dadf3<^{rj), 



'^^^\r,) = d^dp^{r,)^ 



^ *(^)(f,).*(^)(f,). (2) 
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Since V^<i>(7^ = Tr$*^^''(r) — 0, the principal curvatures 
of the pseudopotential cannot be chosen independently. 
Our algorithm requires specifying the electric potential 
curvature tensors ^^'^\rj), which can be determined up 
to their irrelevant sign from the corresponding ^^^-'(rj) 
through Eq. (|2]), provided that is fully confining at rj. 

We can quantify the strength of a microtrap indepen- 
dently of the distance to the electrode plane Zj and the 
potential ?7if cos(rirfi) applied to the rf electrodes by [T2] 



n,^\det^^^\r,)\^/^x{zyu,i). 



(3) 



For typical experimental parameters C/rf — 50 V, firf = 
2tt X 200 MHz, and Zj — 30 /im, the geometric mean of 
the three principal trapping frequencies for ^Be"*" ions is 
u] — Kj X 27r X 53 MHz. The values Kj depend solely on 
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the geometry of the surface electrodes. Small dimension- 
less curvatures Hj can to some extent be compensated 
by lowering fiif, subject to the Mathieu stability require- 
ments |10j . or by increasing U^f. However, much can be 
gained by optimizing the electrode shapes such that the 
Kj are maximized for given constraints. The trap depth, 
i.e., the ponderomotive pseudopotential ^{fg-pj) = Vj'^j 
of the lowest saddle point rgp between a microtrap and 
the boundaries, is not optimized in the present algorithm, 
but will be given in what follows in terms of the energy 
scale = q'^U^f/{4mfl'^fz'j). For the above typical pa- 
rameters, = 4.7 eV. 

The resources we consider for realizing general surface 
electrode arrays for trapped ions are two types of peri- 
odic electrodes in the xy plane, one being grounded or 
at a slowly varying control potential, and the other car- 
rying the rf potential at a uniform amplitude and phase. 
These electrodes are assumed to completely cover a sin- 
gle plane ^TH [TB] . Our algorithm optimizes the shapes 
of these electrodes to achieve maximal curvatures Kj at 
a set of M field-free microtraps per unit cell at posi- 
tions {rj}j'Li and with electric potential curvature ten- 
sors {^^^\rj) = CU,{rj}fLj^. The 3 X 3 matrices 
must be symmetric and traceless to fulfill the Laplace 
equation. The dimensionless curvatures of the various 
microtraps are Kj = \C\ x Zj \ det Tj\^/^. Maximizing the 
common scaling factor |C| therefore simultaneously max- 
imizes all the Kj while preserving the relative curvatures 
of the different microtraps. 

We subdivide the unit cell of the desired Bravais lat- 
tice into N small patches, and label i?i — ^i/Uri the di- 
mensionless electric potential amplitude of the i*^ patch 
electrode, with = (t^^ = 1) signifying a connection 
to ground (rf). The algorithm below produces the op- 
timal set of binary values i?i G {0, 1} Vi = 1 . . . with 
a negligible number of instances where intermediate val- 
ues < -di < 1, which would require additional rf volt- 
age sources, occur. The radio-frequency electric poten- 
tial ^{x,y,0) in the electrode plane is fully determined 
by the shape of the surface electrodes, i.e., by the A^- 
componcnt parameter vector This Dirichlet boundary 
condition can be extended from the electrode plane into 
the third dimension using a Fourier transform [iTJ [13] 
since the electric potential satisfies the Laplace equa- 
tion V^<i>(r) = 0: each Fourier component of in-plane 
wavenumber k — (kx,ky) is damped exponentially as 
exp(ifc2;X -I- ikyy — \\k\\z) away from the electrode plane. 
The typical size of the electrode patches gives a natural 
cutoff in this Fourier series. The resulting dimensionless 
potential Q{§;r) — J2i^i®ii^ then constrained to 
have the desired structure of M microtraps per unit cell: 
at the position of each field-free microtrap the electric 
field must vanish, 

Ve{d; r,) = = 0. (4) 



giving 3M linear conditions on the coefficient vector i). 
Further, the curvature tensors of the electric potential 
at the microtrap positions must match the desired trap 
curvatures, 

i 

(5) 

giving another 5M linear constraints on i9. In princi- 
ple we could add other linear local constraints, such as 
higher order derivatives of the electric field for control- 
ling anharmonicities, to further customize the result to 
the situation at hand; however, doing so might result in 
a decrease of the maximum achievable curvatures. The 
set of 8M linear constraints are denoted A ■■& = Cb in 
matrix representation. We define 7 = A'^ ■ b as their 
unique inhomogeneous solution which is orthogonal to 
the null space of A, computed from the Moore-Penrose 
pseudoinverse A"*". After decomposing i3 = C^ + d' with 
I?' • 7 = 0, the electrode optimization problem is to 

• find the vector -d satisfying A ■ -d' = A-y-^ ■ -d = 
0, with G {0,1} Vz = 1 . . . A^, and maximizing 

|C| = l^^-7l/ll7lP, 

where we have used 7-*- = 1— 77'''/||7|P as the perpendic- 
ular complement of 7. Such integer linear programs are 
typically nondeterministic polynomial-time (NP) hard 
to optimize [2]. Fortunately, relaxing the integer con- 
straints to < "i^i < 1 Vz = 1 . . . yields a linear pro- 
gram [T5] that gives globally optimal results within 0{N) 
time and memory. The reason that relaxed constraints 
are sufficient in practice is that the optimal solution for 
I? contains at most 8M "interior" values, < -di < 1, 
and consequently at least A^ — 8M patch potentials are 
"railed" at either = or "di = 1 and satisfy the original 
binary constraint. This is because the optimal solution 
of a linear program is a "basic" solution; by construction 
all basic solutions fulfill this railing condition [TS]. As 
the resolution of the patch decomposition is increased to 
infinity, the combined area of the fixed number of un- 
railed patches becomes vanishingly small, and the opti- 
mal electrode converges to a solution of the integer linear 
program in the sense that ^{x,y,0)/U,-f S {0,1} for all 
{x,y). In practice, rounding the interior values to or 
1 has very little effect on the properties of the resulting 
potential even for modest grid resolutions. 

We have implemented the algorithm for oblique and 
hexagonal Bravais lattices. The oblique implementation 
uses N — ni X n2 identical parallelogram-shaped patch 
electrodes per unit cell, while the hexagonal implemen- 
tation takes ni = n2 = n and further divides each rhomb 
into two equilateral triangles, giving A^ = patch elec- 
trodes. Optimizations with A^ = 10^ can be done in half 
an hour on a desktop computer. Our results are opti- 
mized with ni = n2 = n = 250, using (2 x 250)^ Fourier 
waves. We have found it unnecessary to use symmetries 
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FIG. 1: Optimized dimensionless curvatures k (black, left 
axis) and trap depths rj (cyan, right axis) as functions of the 
ratio of microtrap height z to inter-ion spacing d, for sev- 
eral lattice geometries. Curvature tensors are cylindrically 
symmetric out-of-plane quadrupoles as for the optimized ring 
trap (k = 0.298 and i) = 0.0196; see text). The decreases 
in dimensionless curvatures and trap depths depend strongly 
on the typical electrode structure size A = 2-k/Q (see text). 
Spurious microtraps may be present for some parameters. 





FIG. 2: Top (bottom) row: optimized electrodes for square 
(triangular) lattices of microtraps at varying height z/d = Q.2, 
0.3, 0.4, 0.5 (0.25, 0.5, 0.75, 1) in units of the ion spacing d. 
Ground and rf electrodes are shown in white and blue (with 
arbitrary assignment). Microtrap locations are marked with 
triangles. Curvature tensors are cylindrically symmetric out- 
of-plane quadrupoles as for the optimized ring trap (see text). 
For z/d = 0.2, 0.3 (0.75, 1) the electrodes maximizing the di- 
mensionless curvature k (outlines shown) give rise to spurious 
trapping sites which have been eliminated by imposing addi- 
tional constraints (see text) reducing k by 0.07% and 0.04% 
(4% and 18%), respectively. 



for specific wallpaper groups in the assignment of the 
patch electrodes, as the solutions found by interior-point 
methods of linear programming do not spontaneously 
break such symmetries, and the gain in computational 
time is only minor, due to the 0{N) scaling. 

The surface-electrode ion trap with the highest possi- 
ble curvature is an rf ring electrode with a ratio of outer 
to inner radii of 4.98 embedded in an infinite grounded 
plane 13J. We use its dimensionless curvature k = 0.298 




FIG. 3: Optimized electrodes for a bilayer honeycomb lattice. 
The desired microtraps, shown as equipotential ellipsoids, are 
at heights z±/h — 0.4 ± i^^, with h the honeycomb lattice 

periodicity and d = the inter-ion distance. Vibrational 

frequencies of the microtraps are {uJx,<Jy,0Jz) = ujo{4>~^ ,1, (/)) 
in the three mutually orthogonal nearest-neighbor directions 
(marked with solid lines and projections; lighter shades in- 
dicate higher frequencies), where (f) — (1 + \/5)/2 is the 
golden ratio. Dimensionless curvatures are K-jk = 0.0022 
and K+/k = 0.020 (differing because of the different heights 
z±). The locations of spurious microtraps are marked with 
small spheres. Thin vertical lines are added for clarity of per- 
spective. 



as a reference value. Regular sparse {z/d < 1/2) lattices 
of such cylindrical microtraps at a given trapping height 
z are most easily constructed from lattices of similar ring 
electrodes, with the inner and outer radii re-optimized for 
each z/d such that k is maximized. This simple param- 
eterization already gives surprisingly good results, with 
values of k at most 10% below the global optima, as 
determined by high-resolution linear programs (Fig. [TJ. 
However, the outer radius must be smaller than d/2, 
which limits the possible trapping heights or necessitates 
more complex parameterizations. Nonparametric opti- 
mization, on the other hand, has solutions for any desired 
trapping height, as shown in Fig.[l] albeit with exponen- 
tially decreasing dimensionless curvatures for large z/d 
[k (X (Qz)^ exp(—Qz) for Qz 3> 1, with Q the wavenum- 
ber of the dominant Fourier component of the electrode 
pattern] . The dimensionless trap depths rj are also found 
to decrease exponentially with z/d [77 oc exp{—2Qz) for 
Qz ^ 1]; however, the given values are not necessarily 
maximal. Figure |2] shows that for low trapping heights 
the electrodes consist of almost circular ring electrodes 
centered under the desired microtrap locations; but at 
larger values oi z/d the electrodes of neighboring lat- 
tice sites overlap and interact significantly. This leads to 
nonintuitive patterns that, by construction, are optimal 
in the sense of maximizing k. 
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A major advantage of our nonparametric method is 
that much more complex trapping geometries can be con- 
structed with no additional effort. However, for compli- 
cated microtrap structures the optimized dimensionless 
curvatures kj and trap depths rjj tend to be vanishingly 
small, requiring unfeasibly large rf potentials for suffi- 
cient trapping. As an example. Fig. [3] shows the op- 
timized electrode for constructing a bilayer honeycomb 
lattice with right angles between nearest-neighbor direc- 
tions [16' and unequal trapping frequencies in these direc- 
tions. The microtrap frequencies are chosen to be equal 
in the two trapping planes, leading to different dimen- 
sionless curvatures k±. 

Figures |2] and |3] illustrate that optimal electrodes tend 
to consist of large smooth structures, with the un-railed 
patches located on the boundary between the rf and 
ground electrodes, reminiscent of anti-aliasing effects. 
This observation, which ensures that optimal electrodes 
can realistically be microfabricated, is attributed to the 
fact that microtraps can often be approximated by a lo- 
calized cancellation of a small number of long- wavelength 
Fourier modes; the amplitudes of long-wavelength modes 
are maximized when the electrodes consist of large uni- 
form areas, with sizes related to the wavelengths of the 
modes involved in the microtraps. 

The linear programming algorithm is not limited to in- 
finite periodic planar gapless electrodes. Any set of two- 
or three-dimensional electrodes can take the role of the 
"patch" electrodes, with an appropriate generalization of 
the method for extending the electric potential into the 
region where the ions are to be trapped. The algorithm 
will then specify which of these electrodes should be con- 
nected to rf and which to ground. In particular, finite 
planar electrodes for a typical experimental setup may 
be optimized in order to avoid finite-size effects coming 
from the truncation of a periodic electrode array. 

The algorithm offers unlimited freedom for placing ar- 
bitrary microtraps in arbitrary locations, and is optimal 
in the sense of globally maximizing the curvatures ([3| 
and producing smooth binary electrode shapes, but also 
has some limitations. First, it offers no possibility to 
optimize the trapping depth instead of its curvature, be- 
cause the former is a spatially nonlocal property and is 
therefore incompatible with the linear programming al- 
gorithm. Possibly parametric optimizations such as the 
ring-trap optimizations outlined above can be adapted 
to take the trap depth into account [HI [13]. Second, 
specifying the desired set of microtraps does not imply 
that these will necessarily be the only microtraps gener- 
ated by the optimized electrode set. Complex electrodes, 
such as those shown in Fig. [3] often generate a multi- 
tude of unwanted microtraps, predominantly on points or 
lines of symmetry of the underlying wallpaper symmetry 
group. In highly symmetric setups (Fig. [2]), these spuri- 
ous microtraps can often be eliminated by constraining 
the out-of-plane electric field Ezir) to judiciously chosen 



nonzero values at the spurious trap sites, thus adding fur- 
ther linear equalities or inequalities to the algorithm and 
reducing the optimized curvatures Kj. In setups lacking 
symmetry (Fig. [3]) this procedure is much less intuitive. 
In many cases it is possible to operate the trap in a regime 
where the spurious trapping sites are unstable, while the 
trapping sites at rj are stable [TU] . 

The Coulomb interaction between ions at different 
trapping sites can lead to net static forces, due either 
to a finite lattice or to an insufficiently symmetric or 
non- planar unit cell (giving out-of-plane forces in Fig.|3|. 
Such static forces tend to push the ions away from their 
rf zeros, thus inducing unwanted micromotion. Efforts 
should be made to compensate these forces by suitable 
control potentials. For this purpose we can subdivide the 
unit cell into control electrodes and "patch" electrodes 
before optimizing the latter with our algorithm; alter- 
natively, the optimized ground and rf electrodes can be 
subdivided into areas biased with different control volt- 
ages. The same approaches are readily generalized to 
non-periodic ion trapping setups. 
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